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In a two-dimensional quantum dot in a GaAs heterostructure, the spin-orbit scattering rate is 
substantially reduced below the rate in a bulk two-dimensional electron gas [B.I. Halperin et al, 
Phys. Rev. Lett. 86, 2106 (2001)]. Such a reduction can be undone if the spin-orbit coupling 
parameters acquire a spatial dependence, which can be achieved, e.g., by a metal gate covering 
only a part of the quantum dot. We calculate the effect of such spatially non-uniform spin-orbit 
scattering on the weak localization correction and the universal conductance fluctuations of a chaotic 
quantum dot coupled to electron reservoirs by ballistic point contacts, in the presence of a magnetic 
field parallel to the plane of the quantum dot. 
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In disordered metals, interference of time-reversed tra- 
jectories leads to a small negative correction to the con- 
ductivity, known as weak localization.Ell3 With strong 
spin-orbit scattering, the effect of such interference is op- 
posite, causing weak aniilpcalization, a positive correc- 
tion to the conductivity&Q Both interference corrections 
are suppressed when time-reversal symmetry is broken by 
a magnetic field. The same phenomena are observed in a 
two-dimensional electron gas, such as is formed in GaAs 
heterostructures. While spin-orbit scattering in metals 
is largely due to scattering from the metal ions or from 
impurities, in a GaAs heterostructure, spin-orbit effects 
mainly arise from the asymmetry of the potential creat- 
ing the quantum well (Rashba term), as well as from the 
lack of inversion symmetry which may occur in the crys- 
tal structure of the material forming the heterostructure 
(Dresselhaus term). 

Recently, it has become possible to study spin-orbit 
scattering in finite size systems, -such as metal grains and 
semiconductor quantum dots.ETEl In the universal regime, 
where all relevant time scales (spin-orbit time r so , inverse 
level spacing/broadening) are much larger than the elec- 
tron transit time T crg , such|-systems can be described us- 
ing random-matrix theorya Even though r so 3> r org in 
the universal regime, spin-orbit scattering may still have 
a significant effect on wavefunctions and transport prop- 
erties if t so is comparable to the inverse level spacing or 
level broadening, respectively. 

For metal grains, the spin-orbit Hamiltonian H so is 
modeled tw-,a random hermitian matrix with symplectic 
symmetryji!] the same symmetry as in the case of bulk 
disordered metals. However, for GaAs quantum dots, 
the situation is more complicated: both the symmetry of 
the random matrix representing _ff so and the spin-orbit 
time are diflepept from the case of a bulk two-dimensional 
electron gasJlalia The complications arise from the spe- 
cial form of the spin-orbit Hamiltonian H so in GaAs, 
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where A^ and X y are length scales describing the-spin- 
orbit scattering strength in a GaAs heterostructurelI3 and 
a are the Pauli matrices. The structure of H so is that of a 
"non-Abelian vector potential" , coupled to the electron's 
spin.Ed As this "vector potential" has no spatial depen- 
dence, and hence no "flux" , a suitable gauge transfor- 
mation removes the spin-orbit scattering term from the 
Hamiltonian up to corrections of order LjX XiV which arise 
due to the non-Abelian nature of the "vector potential" 
of Eq. ([!]). (Here L is the size .of the dot; In the univer- 
sal regime one has L <C A x . a .E£l) As a result, the spip. 
orbit scattering time is increased by a large factoiErEJ 
~ X x X y /L 2 ~ T^/r CTg over its value t££ ~ 2\ x X y /£vF in 
a bulk two-dimensional electron gas with Fermi velocity 
vf and mean free path £ equal to that in the dot, or with 
i ~ L for the case of a ballistic dot. Moreover, as was 
shown in Ref . the symmetry of the transformed spin- 
orbit scattering term is not symplectic, but unitary.E3 

In this paper we investigate the case where the spin- 
orbit coupling parameters A^ and X y are not constant 
throughout the quantum dot. Experimentally, such a 
situation could be created with the help of a metal gate 
parallel to the two-dimensional electron-ass that changes 
the asymmetry of the quantum well.tSO If the metal 
gate covers only a part of the quantum dot, as is shown 
schematically in Fig. El the translational invariance of 
H so is lifted. Hence, the spin-orbit scattering can no 
longer be gauged away to leading order in L/X x , L/X y . 
In other words, the "non-Abelian vector potential" in Eq. 
(^) now represents a nonzero "flux" . The consequence is 
a significant increase of the spin-orbit scattering rate and 
a restoration of the symplectic symmetry of H so . Thus, a 
metal gate that changes the asymmetry of the quantum 
well in only part of the quantum dot has a fundamen- 
tally different effect than a metal gate that changes the 
quantum well potential uniformly throughout the dot. A 
gate that covers only part of the dot is clearly the more 
effective tool to tune the spin-orbit scattering rate. 

As an example, let us consider a quantum dot for 
which the spin-orbit scattering is due to the asymme- 
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try of the potential well (Rashba term) only, so that 
I A x | = | A j, | = A.cj If there is a metal gate over half of the 
dot as in Fig. EL the spin-orbit scattering rate may take 
two different values A -1 = A -1 ± ^X~ x in the two halves 
of the dot. As discussed in Refs. fllfL2l the spatially uni- 
form component of H so leads to a unitary perturbation 
of the Hamiltonian, with a characteristic time 



'erg 



(A/L) 4 ~ (C) 2 /Te: 



The spatially varying component gives rises to a symplcc- 
tic perturbation of the Hamiltonian, with a characteristic 
time that can be estimated as the time to accumulate a 
"flux quantum" from the "vector potential" in Eq. (flu , 
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If A s ~ A, t s s becomes comparable to t££ , the spin-orbit 
scattering time in a bulk two-dimensional electron gas. 
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FIG. [j]. Schematic drawing of a quantum dot with a metal 
gate (hatched) over part of the dot. The role of the gate is to 
change the asymmetry of the potential of the quantum well 
beneath it, and hence the spin-orbit parameters \ x and \ y . 

We now present a quantitative calculation of how such 
a tunable spin-orbit scattering time affects the quantum 
interference corrections to the conductance: the weak 
(anti) localization correction (SG) and the conductance 
autocorrelation function cov[G(B),G(B')], where B is 
a magnetic field. Qitr work extends previous works of 
EfetovEj and FrahmEj for the magnetic-field dependent 
quantum interference corrections without spin-orbit scat- 
tering or a parallel magnetic field. The effect of the spa- 
tially uniform component of H so on the weak localization 
correction was calculated in Ref. [l^. 

We consider a quantum dot coupled to two electron 
reservoirs (labeled 1 and 2), via ballistic point con- 
tacts that have N% and N 2 channels each. We assume 
that the electron motion in the quantum dot is chaotic, 
so that random matrix theory can be used to calcu- 
late the conductance distribution in the universal regime 
gHBB,Ti/T so <C h/r org .U (Here g is the electron g-factor 
and hb the Bohr magneton.) The quantum dot is de- 
scribed in terms of its scattering matrix S, which, for 
particles with spin, is a N x N unitary matrix of quater- 
nions, N = Ni+N 2 . Quaternions are 2x2 matrices with 



special rules for transposition and complex conjugate. 
Starting point of the calculation is the Landauer formula 
for the two-terminal conductance G of the quantum dot 
at zero temperature, 



G= 2 4^-^rSAS^A. 
h N h 
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(2) where the diagonal matrix A has elements 



A, 



N 2 /N, j = l,...,N u 
—Ni/N, j=N 1 + l,.. 



In order to find the average and variance of the con- 
ductance G it is sufficient to compute the average 

(S khtlu {s,B)Sm<^ v >{s',B'Y) 

in the presence of spin-orbit scattering and for arbitrary 
values of the magnetic field B and Fermi energy e. (Ro- 
man indices refer to the propagating channels in the 
leads, greek indices refer to spin.) In a random-matrix 
approach, the statistical properties of the scattering ma- 
trix S can either be calculated from a hermitian random 
matrix that represents the Hamiltoniaa of the quantum 
dot, or from a random unitary matrixB Here we use the 
latter approach; equivalence of the two approaches, in- 
cluding the dependence on an external parameter, was 
shown in Ref. Ell The N x N matrix S is written aaEJ 
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where U is an M x M random unitary symmetric-matrix 
taken from Dyson's circular orthogonal ensembleo and R 
is a unitary matrix of size M — N. The N x M matrix P 
and the (M — N) x M matrix Q are projection matrices 



with Pij = Si j and Qi 
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The quaternion ele- 



ments of the matrices U, P, and Q are all proportional 
to the 2x2 unit matrix 1. The matrix R is given by 



R(s, B, r so ) = exp 
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where A is the mean level spacing of the dot and H' is an 
[M — N) dimensional quaternion matrix generating the 
perturbations to the dot Hamiltonian that correspond to 
the magnetic field and the spin-orbit scattering, 
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Here Aj (j = 1,2) and X are real antisymmetric ma- 
trices of dimension M — N, with tr AiAj — M 2 Sij and 
tr X X T — M 2 .El The symmetry of the spin-orbit term in 
Eq. (^) is chosen in accordance with Eq. ([!]) , taking into 
account that the spin-orbit Hamiltonian has symplectic 
symmetry once the coupling parameters \ x and X y de- 
pend on position. In Eq. (In) the orbital and Zeeman 
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effects of the magnetic field have been separated. The 
first term describes the Zeeman coupling to the spin of 
the electrons. The second term models the orbital effect, 
where x is related to the perpendicular component of the 
magnetic field, 

x 2 = ce 2 L 4 Bl/{hT crg A). 

c being a numerical coefficient of order unity! 2 ! At the end 
of the calculation, the limit M — > oo should be taken. 

We now describe our calculation, which was done to 
leading order in l/N. Corrections for finite A are dis- 
cussed at the end of this paper. To leading order in 1/M 
and l/N, it is sufficient to consider the elements of U 
as random Gaussian variables with zero mean and with 
variance (UijUjg) = M~ 1 (5 i kS : ji + SuSjk)£3 We then ex- 
pand Eq. (H) in powers of U and perform the Gaussian 
averages to leading order in 1/M. We thus find 
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where, in tensor notation, 

D = (Ml (g> 1 - tr R ® i?' t ) _1 , 

C = {Ml®t-tvR®R'*y 1 . 

Here R'* is the quaternion complex conjugate of R', 
T = 1 !g> <T2) and the tensor multiplication should be 
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(G(e,x,b)) = 



2e 2 AiA 2 
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understood as "backwards multiplication" for the second 
matrix, i.e., with the multiplication rules 

((Tj <g) CTj)(0V ® CTf) = {(Ji(Ji>) ® (crj'(Tj). 

The two contributions C and D are the equivalents of 
cooperon and diffuson in the conventional diagrammatic 
perturbation theoryo 

Taking the limit M — > oo and defining a dimensionless 
magnetic field b — ngfisB / A and spin-orbit scattering 
rate a 2 = 27t?i/t so A, we find 



D^ 1 = N D (t ® 1) +ib- ((7 1)- ib' ■ (1®ct) 
+ 2a 2 (l (g> 1) - a 2 (cr 2; ® (Tj; + cr y (g> CT y ), 

G -1 = A c (l (8>l) + ib- (a® 1) + ib' • (1 ® a) 
+ 2a 2 (l ® 1) - a 2 (cr x ® + cr,, ® cr y ), 

where A^ and Ac are given by 

Ad — N — 27ri(e - e')/A + (1/2) (a: 
Ac = N — 2ni(e - e')/A + (1/2) (sc - 
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We now set b = bx, b' = b'x, take the inverses in Eqs. 



(|9|) and (10), and calculate the average and covariance of 
the conductance G from Eq. (|). For (G), we find 
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A c + 2a 2 A c + 4a 2 N C (N C + 2a 2 ) + 4b 2 



(13) 



where Ac = AT + 2x 2 , as follows from Eq. ( [i"2"| ) with To calculate the zero temperature conductance 

fluctuations, it is sufficient to know the two-point correlator Ja) to leading order in l/N. (Contributions from higher- 
order correlators vanish since they contain a factor tr A = OtJ) We then find 
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2a 4 + 2(b + b') 2 + 2|A D + 3a 2 | 2 

| (b - b') 2 + A D (2a 2 + N D )\ 2 T + b') 2 + {N D + 4a 2 ){N D + 2a 2 ) | 



Fc is obtained from Fp by the substitution x 1 — > —x', 
b' — ► — b', and Ad — > Ac, and Ad and Ac are given in 
Eqs. ( |li"| ) and ([l2|) above. The conductance fluctuations 
at finite temperature are obtained from Eq. ( |l4| ) by mul- 
tiplication with the derivatives of the Fermi function at 
energies e and s' and subsequent integration over e and e'. 
Equations (|l3|) and ( |l4] ) recover the result of Refs. [T^Jl9| 
for the spinless case a = b — 0. Further, one has (SG) — 
(e 2 /h)(N 1 N 2 /N 2 ), varG = 2(e 2 //i) 2 (Ai A 2 /A 2 ) 2 if there 
is no parallel magnetic field, while spin orbit scattering 
is strong (a 2 3> A), in agreement i^ith known results for 
the circular symplectic ensemble.u (For comparison, no 
positive (SG) is observed for the case of a spatially uni- 
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form spin-orbit coupling, see Ref. [L2|.) Equations fll^ , 
and (|lj) are illustrated in Fig. || where we show (SG) 
and varG as a function of the perpendicular magnetic 
field x for various values of the dimensionless spin-orbit 
scattering rate a. 

In the presence of both strong spin-orbit scattering 
(a 2 > A) and a large parallel field (b 2 > A 2 , Aa 2 ), 
all terms contributing to the cooperon G are suppressed. 
As a result, there,is no weak localization correction to the 
conductance,c3'E3 (SG) = 0, and the conductance fluctu- 
ations are reduced by an additional factor two, var G = 
(e 2 //i) 2 (AiA 2 /A 2 ) 2 , see also Ref. [□]. Furthermore, the 
conductance is no longer symmetric under reversal of the 
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perpendicular magnetic field. This observation is best 
illustrated by the correlator {(G(x) — G(—x)) 2 ), which 
goes to 2 var G for perpendicular magnetic field strengths 
x 2 ^> N, see Eq. (|l4|). For comparison, in the presence 
of either spin-orbit scattering or a parallel field (but not 
both), one has G(x) = G(—x) for all x. 

On a phenomenological level, dephasing can be added 
to the cuttfpt description via the voltage probe model of 
ButtikerEa-LJ In this model, a fictitious voltage probe is 
attached to the quantum dot. Electrons escape from the 
dot into the voltage probe at a rate l/r^, where is the 
dephasing time, and are then reinjected from the voltage 
probe without phase memory. The escape into the volt- 
age probe is described by an imaginary term i?i/2T0 in 
the Hamiltonian. With the form (Q) of the Landauer for- 
mula, the reinjection of particles from the voltage probe 
has no effect on the conductance to order N . Hence, to 
leading order in l/N, inclusion of dephasing amounts to 
the replacement 
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FIG. 0. Left panel: Weak localization correction (SG) as 
a function of the dimensionless perpendicular magnetic field 
x, for a — (bottom curve), a = 0.3, a — 0.5, a = 1.0, and 
a — > oo (top curve). Right panel: var G versus x for the same 
values of a (a — is top curve, a — > oo is bottom curve). In 
both cases, there is no parallel component of the magnetic 
field and G is measured in units of (e 2 /h)N\N2/N 2 . 
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